Dissipative defect modes in periodic structures 
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We show that periodic optical lattices imprinted in cubic nonlinear media with strong two-photon 
absorption and localized linear gain landscapes support stable dissipative defect modes in both fo- 
cusing and defocusing media. Their shapes and transverse extent are determined by the propagation 
constant that belongs to a gap of the lattice spectrum which, in turn, is determined by the relation 
between gain and losses. One-hump and two-hump dissipative defect modes are obtained. 



PACS numbers: 

Implementation of stable spatial localization of light 
is a topic of continuously renewed interest because of 
its fundamental importance and potential applications. 
The use of structures with shallow periodic refractive in- 
dex modulations is a fruitful strategy to generate and 
to control soliton-like light packets (see, e.g. 0, for 
recent reviews). Localization, alias defect modes, can 
be achieved in linear periodic systems with structure de- 
fects. The interplay between nonlinearity, periodicity of 
the underlying structure, and local lattice deformations 
results in appearance of new types of localized modes and 
important modification of stability of modes existing in 
linear regime (see [3|-[6| in the context of nonlinear optics 
and in [JQ in the context of Bose-Einstein condensates). 
Further development of this concept can be achieved in 
non-conservative systems where dissipation and gain are 
present, in particular, using dissipative defects. The cor- 
responding idea, where the defect is produced by local- 
ized gain, has been suggested upon study of interactions 
of gap solitons in shallow fiber Bragg gratings [9| and in 
media of two- level atoms , as well as for stabilization 
of spatial dissipative solitons described by the complex 
Ginzburg-Landau equation [11] (i.e. relevant to cavity 
solitons). Our statement is also related to the studies of 
the discrete Ginzburg-Landau equation, which possesses 
unstable localized solutions [ijj], one of the stabilization 
mechanisms for which is the quintic nonlinearity (l3| . 

In this Letter we report existence and stability of dissi- 
pative defect modes in periodic optical lattices (we distin- 
guish them from solitons since the modes are pinned by 
several " defect" channels where gain is realized and exist 
only due to competition between the gain and nonlinear 
losses). Some of the properties of such modes appear to 
be very different from their conservative counterparts 0- 
,8]. Being attractors they possess enhanced stability and 
can be excited from sufficiently large class of initial con- 
ditions. 

Specifically, we consider propagation of laser radiation 
in a lattice with spatially localized linear gain and nonlin- 
ear losses (two-photon absorption) that can be described 
by the nonlinear Schrodinger equation for the dimension- 



less light field amplitude q: 
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q vv - R{rj)q + ij(v)q 



a\q\ 2 q-ia\q\ 2 q (1) 



Here 77 and £ are the normalized transverse and longi- 
tudinal coordinates, respectively, a = 1 (a = —1) for 
focusing (defocussing) media, a > characterizes nonlin- 
ear losses, R(rj) is a periodic function, and 7(17) describes 
the linear gain which is nonzero only in a finite domain, 
i.e. 7(77) = for \rj\ > 770 (770 is a positive constant). 

One of our central results is that system ([T]) supports 
stable localized modes in both focusing and defocusing 
nonlinear media. The complex amplitude of such modes 
can be written in the form 77) = e lb ^ +ld ^w{vi), with 
a real propagation constant 6, real stationary amplitude 
w and phase 9, and with the zero asymptotes: w(rj) — > 
at 77 — > ±00. Consequently we rewrite (JXJ) in the form 
of equations for real- valued functions 10(77) and j{rf) = 
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hw= ^Tl_ J_ + aw A + Rw ^ ^ = 2 1W 2 - 2aw i (2) 



The propagation constant of an exponentially localized 
mode must belong to a gap in the spectrum of the peri- 
odic "potential" R(rj), i.e. to an empty domain in Fig. [I] 
(a). This is another important finding of the present 
work, which suggest a way to classify defect modes. In- 
deed, let us consider the limit 77 — > 00. From the second 
of Eqs. @ it follows that j = 2a f°° w 4 dr] for 77 > 7j . 
Hence, (i) there exists an energy flow outwards the impu- 
rity, j > 0, and (ii) the following relation (obtained us- 
ing l'Hopital's rule) lim, ) _ ! . 00 j/w 2 = — a lirn^oo w /w v 
holds. For localized modes with the finite energy flow 
U = J_ oo w 2 (rj)dr], i.e. with w 2 decaying more rapidly 
than I/77, one obtains that the above mentioned limit 
is zero as long as w v is sign-definite. Otherwise, if w v 
changes sign at large 77, one also verifies that the estab- 
lished condition j > can hold only if the above men- 
tioned limit is zero. Hence, j 2 /w 3 decays faster than w 
and can be neglected at rj — > 00 and b satisfies the Hill 



equation bw 



R(r))w. Thus, to guarantee the 
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FIG. 1: (a) The band spectrum of R(r)). For p r = 5 the lower 
edge of the semi-infinite gap is at b — 2.875 while the first 
finite gap corresponds to the interval (—0.645, f.840). Shad- 
owed (empty) regions indicate the allowed bands (gaps). De- 
pendence U vs pi for the one-hump (curve I) and two-hump 
(curve 2) modes is shown in (b) for the focusing and in (c) 
for defocusing media, (d) Propagation constant vs gain for 
one-hump modes in focusing (curve I) and defocusing (curve 
2) media. Circles in panels (b)-(d) correspond to the modes 
shown in Fig. [2] In all the cases a = 1. 



realized only in one period, then 7(17) = Pi cos 2 (Sir)) for 
\v\ < ?7o = 7r/(2Q) and 7 = for \r)\ > t]q. Notice, that 
other functional profiles of the gain in amplifying lattice 
channels yield qualitatively similar results. Further we 
set SI — 2, p r — 5 and vary pi and a. 

When the propagation constant belongs to the semi- 
infinite gap, the defect modes can be obtained in focusing 
medium only. The emerging branches of the modes are 
shown in Fig. Q] (b), while examples of their profiles are 
presented in the left column of Fig. [5J As we mentioned, 
b is determined by the gain/loss as shown in Fig. [TJd). 
Notice that depending on whether the gain is applied to 
one or two channels the attractor is a one-hump [panels 
2(a), 2(b)] or a two-hump [panel 2(c)] mode. 




decaying asymptotic of the mode, b must belong to a gap 
of the lattice spectrum. We emphasize the importance of 
the sufficiently fast decay of the linear dissipation coeffi- 
cient for the last conclusion. 

In contrast to conservative defect modes, the station- 
ary states in our system exist not only due to balance 
between nonlincarity and refraction/diffraction of the lat- 
tice, but also due to balance between the linear gain 
and nonlinear losses ^JUl, which follows from Eq. ((T|): 
/- l{i))w 2 dr) = a J_ w drj. Therefore, the propaga- 
tion constant b and the energy flow U are determined by 
linear gain 7 and by the nonlinear losses a. 

Turning to numerical study of the modes, we chose 
R(rf) — p r cos 2 (Sir)) to model the lattice, whose spec- 
trum is shown in Fig. [IJa). We consider localized linear 
gain profiles, i.e., we assume that 7(77) oc R(r)) on one or 
several periods of R(t)). For instance, if amplification is 



FIG. 2: Profiles of one- hump (one-channel gain) (a),(b), (d), 
(e) and two- hump (two-channels gain) (c), (f) modes in focus- 
ing (left column) and defocusing (right column) media. The 
humps in two-hump modes in focusing and defocusing media 
are out-of-phase and in-phase respectively. 

For b in the first finite gap we have found the stable 
defect modes in the defocusing medium. Their branches 
are shown in Fig. Q] (c), and the examples of the ampli- 
tude configurations are given in the right-column panels 

of Fig. m 

Like in the conservative case, when the propagation 
constant approaches a gap edge defect modes extend far 
beyond the lattice channels where amplification is ap- 
plied, exhibiting pronounced shape oscillations which ap- 
proaches properly normalized Bloch states w bordering 
the respective gap edge. At fixed a in a focusing medium 
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an increase of the gain pi results in progressive increase 
of the mode peak intensity and in gradual contraction 
of light into the amplifying channel [Fig. Ufa), (b)]. In 
defocusing medium the peak amplitude also grows with 
Pi but the modes possess strongest localization at in- 
termediate values of pt and may extend dramatically 
across the lattice when the gain intensity becomes too 
high [Fig. |2fd),(e)]. The energy flow U monotonically 
grows with pi [Fig. tUb), ( c )]- For sufficiently small and 
moderate nonlinear losses one-hump modes exist even at 
Pi — > in both focusing and defocusing media. In this 
limit b approaches the band-edge, the amplitude of the 
mode becomes small and its width grows dramatically. 
The nonlinear losses diminish as well, and the balance 
between gain and losses requires the linear gain to go to 
zero. Notice, however, that for sufficiently high nonlin- 
ear losses stable defect modes can be found only above 
certain critical value of pi and in this parameter range 
stable modes exist above threshold energy flow. 

Higher-order modes possess energy flow threshold for 
their existence [curves 2 in Fig. [ljb), (c); we confirmed 
this for all modes in focusing (defocusing) medium that 
change (do not change) their sign between neighboring 
lattice sites] . Increase of pi results in monotonic growth 
of b in focusing medium accompanying mode contrac- 
tion [Fig. Did), curve 1] and decrease of b in defocusing 
medium resulting in mode expansion across the lattice 
when b reaches the lower edge of the gap at sufficiently 
large k [Fig. [TJd), curve 2]. 

Representative examples of the field modulus distri- 
bution of defect modes are shown in Fig. [5J As it was 
mentioned above, when amplification is realized only in 
one lattice channel one obtains simplest one-hump modes 
[Figs. |2{a),(b)]. When two central channels are ampli- 
fying we found two-hump modes in focusing [Fig. [2fc)] 
and in defocusing [Fig. Eff)] media. These states were 
obtained by direct simulation of Eq. (fTJ) with different 
inputs (ranging from localized input beams to extended 
noisy distributions) up to a huge distance. This guar- 
antees that final states that do not change upon propa- 
gation are stationary attractors of Eq. (UJ [solutions of 
Eq. (fT]) with other symmetries are possible too, but they 
are unstable and not shown here] . Notice that for larger 
number of amplifying channels more complicated modes 
with larger number of humps can be obtained. 

The domains of existence of dissipative lattice modes 
on the plane (a,pi) are shown in Fig. [3l In the focusing 
medium one-hump modes exist for all values of pi above 
curve 1 in Fig. Efc) while two-hump modes are found for 
Pi values above curve 2 (the minimal value of pi necessary 



for existence of two-hump modes grows almost linearly 
with a). The situation is more complicated in defocus- 
ing medium where one-hump [Fig. EJa)] and two-hump 
[Fig. [21b)] modes exist only in a finite domain of pi be- 
tween curves 1 and 2. The presence of upper limit for the 
linear gain is connected with derealization of high-power 
modes whose propagation constants approach lower edge 
of first finite gap with increase of pi. The domain of ex- 
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FIG. 3: Domains of existence for stable one-hump and two- 
hump modes (boundaries denoted by circles and pentagons, 
respectively) in defocusing (a), (b) and focusing (c) media. 
In the defocusing medium either one-hump (a) or two-hump 
(b) modes modes exist between the curves 1 and 2. In the 
focusing medium one-hump and two-hump modes exist above 
the curves 1 and 2, respectively. 



istence of two-hump modes in terms of pi values expands 
monotonically with increase of a [Fig. Gib)]. 

Summarizing, we showed that stable defect modes can 
exist in optical lattices imprinted in nonlinear medium 
with strong two-photon absorption provided that gain is 
realized in one or several lattice channels. 
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